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A proposal for applying nonadiabatic geometric phases to quantum computing, called double-loop 
method [S.-L. Zhu and Z. D. Wang, Phys. Rev. A 67, 022319 (2003)], is demonstrated in a liquid 
state nuclear magnetic resonance quantum computer. Using a spin-echo technique, the original 
method is modified so that quantum gates are implemented in a standard high-precision nuclear 
magnetic resonance system for chemical analysis. We show that a dynamical phase is successfully 
eliminated and a one-qubit quantum gate is realized, although the gate fidelity is not high. 
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I. INTRODUCTION 

Precise and reliable control of a quantum system is an 
attractive and challenging experimental issue in quan- 
tum physics. In particular, the importance of its appli- 
cation to quantum information processing has been in- 
creasing recently. A promising way to achieve this is 
to employ geometric phases (or, more generally, non- 
Abelian holonomies) [l|, [2| , because geometric phases are 
expected to be robust against noise and decoherence 
under a proper condition [3|, [4( . A large number of 
studies for applying their potential robustness to quan- 
tum computing have been done, e,g., phase-shift gates 
with Berry phases [5|, nonadiabatic geometric quantum 
gates [6 J? 
puting 14 



, fiol, ITlL [l2l, [l3| , h olonomic quantum com- 
lE EE E3, El, Ella, El| , quantum gates with 



noncyclic geometric phases [22|, and so on. 

For physical realization of geometric quantum gates, 
elimination of a dynamical phase is a key point. Jones 
et al. [5] implemented a controlled phase-shift gate with 
a Berry phase [23|, l2J] by quasi-statistically, or slowly, 
controlling an effective field in a rotating frame. They 
nulled dynamical phase effects using a conventional spin 
echo approach [25]. Zhu and Wang [8] proposed a geo- 
metric quantum gate based on the Aharonov-Anandan 
phase [2y, [gTj , which should be fast since a nonadiabatic 
process is employed. In their proposal, elimination of a 
dynamical phase is achieved by a double-loop method, 
where a dynamical phase cancels out while a geometric 
phase accumulates along two loops. 

Although several experimental techniques for the ap- 
plication of geom etr ic p hases to quantum computation 
are available]!]! E3, EM 5 explicit implementations of ge- 
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ometric phase gates have not been extensively studied 
so far. Without explicit implementations, the often-cited 
advantage of the holonomic quantum gates is nothing 
more than a desk plan. Thus, such explicit examples 
are highly desirable. In this paper, we combine Zhu 
and Wang's approach with Jones et aZ.'s one, employ- 
ing an Aharonov-Anandan phase for fast gate operation 
and a spin echo technique for dynamical phase cancella- 
tion, and demonstrate one-qubit gates with a commercial 
liquid-state nuclear magnetic resonance (NMR) system. 
In many experiments of nonadiabatic geometric quan- 
tum gates [Hi, ll2|, ll3[ , the gate operations in which the 
dynamical phase is arranged to vanish 0, HH have been 
adopted. In the present paper, we show that we may 
have another option for physical realization of geometric 
quantum gates. 

The next section summarizes Zhu and Wang's theoret- 
ical proposal and our modifications for practical imple- 
mentation in liquid-state NMR quantum computer. In 
Sec. lIIII we describe details of our experiments, where im- 
plemented gates are evaluated by performing quantum 
process tomography. Sec.[lVlis devoted to summary. 



II. THEORY 

A. Quantum gates with orthonormal cyclic vectors 

The Aharonov-Anandan phase is a geometric phase 
associated with nonadiabatic cyclic time evolution of 
a quantum system [H, (2?J. Let us write a state at t 
(0 < t < t) as \ijj(t)) in the Hilbert space TL with di- 
mension n. The time evolution of a system is given by 
the Schrodinger equation with a time-dependent Hamil- 
tonian H(t). We take the natural unit in which h = 1. 
The nonadiabatic cyclic time evolution of the quantum 
system means that 

\ 1 p(r)) = U(r)\m)=e i ' t \m), 
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where U(r) is the time evolution operator and 7 G M. 
Let us call 1^(0)) a cyclic vector [26j. We note that the 
dynamical phase 7d associated with this time evolution 
is 

7d = - f \m\H{t)\m)dt, (i) 
jo 

while the geometric component is defined as 
7g =7-7d- 

We take a set of orthonormal cyclic vectors, {|^/c(0))} so 
that 



K( T )> 



u*(0)>, 



where k — l,...,n. Hereafter, we write 1^(0)) as \u k ) 
for brevity. A general state |^(0)) G TL is expressed 
as |^(0)) = Y2=i a k\ u k)- Then, we have |*(r)) = 
E/(r)|tf(0)) = YT k =iake ilk \u k ). We denote a fixed ba- 
sis in TL as {|Z)}, which corresponds to the computa- 
tional basis {|0), |1)} for the case n = 2. In terms of 
{|Z)}, we have |^(r)) = ^^^Tq 1 c^(t-)|Z). It means that 

Q>k = (ufc|*(0)) = Z)r=o lc «( u fel0j where q = q(0). 
Therefore, we find that 



n — 1 n 



n-1 



c«(r) = e^(ZK>( Ufe |0^ = J]) VW, (2) 

r=ofe=i r=o 



where 



=53e < T*(ZK>(« fc |Z , >. 



(3) 



/C=l 



We have j k = 7/c,d + 7/c, g , where 7^ and j k , g are the 
dynamical and the geometric phases associated with \uk), 
respectively. 



B. Single-qubit case 

Let us consider a single-qubit system. We choose the 
Hamiltonian H\(t) as 



Hi(t) 



with NMR in mind, where a = (a x , <j y , <j z ) and 



Q(t) = (coi cos(cu r ft 



-lui sin(cc; r ft - 0), cjq) • (4) 



We note that ft(t) is a time dependent parameter cor- 
responding to the external field and a a (a = x, y, z) 
is the ath component of the Pauli matrices. One can 
freely control c^o by taking a proper rotating frame. 
The transformation to the rotating frame with the fre- 
quency o; r f is made possible by the unitary transforma- 
tion U — e lUriazt / 2 and the Hamiltonian in the rotating 
frame is 



H lT = U ] H X U 



-Qm ■ a. 



(5) 




FIG. 1: Example of dynamics of a single-qubit cyclic vector, 
(a) A time-dependent external field Q(t) and (b) a closed tra- 
jectory on the Bloch sphere corresponding to a cyclic vector 
|-0 + (t)), < t < r = 2ir/\u; T {\. The end point of each arrow 
represents the initial value. We set the parameters uo = 2tt, 
UJl = 0.5 x 27T, cj rf = 0.8 x 27T, and = in Eq. (gj. 

where Q = (u\ + A 2 ) 1/2 , A = ujq - cj rf , m 
(sin \ cos o. sin x sin 0, cosx), and tanx = 0J\/ A. The 
solution of the Schodinger equation is 

\i/;(t)) = e™ rfto "/ 2 e^ m '* /2 M0)). (6) 

We denote the eigenstates of m • a with eigenvalues ±1 
as \ip±). Their explicit forms are 



IV>+> 

iv>-> 



e -**/2 COs||0> 



■e^ /2 sin 111) 



^/ 2 sin^|O)+e i0/2 



cos||l), 



where |0) and |1) are the eigenstates of a z with eigen- 
values +1 and —1, respectively, and taken as the compu- 
tational basis. The corresponding Bloch vectors of \ijj±) 
are 

(V>±|o#±) ±m - 

We require that \ip±) be cyclic vectors. Since 
and are mutually orthogonal, they are identified as 
in Sec. Ill A[ It follows from Eq. ([6]) that the ex- 
ecution time r must satisfy the condition |o; r f|r/2 = 7r, 
i.e., r = 2ir/\uj r f\. Then, \ip±(r)) is written as 



IV>±(t)> 



where j± = tt ± 7r£l/\w T f |. Figure [T] shows an example 
of ft(t) and the closed trajectory on the Bloch sphere 
corresponding to |^+)- Using Eq. (pQ), we find that the 
dynamical phase is 



7±,d = ±-(a;i sinx 



cj cos x) 
while the geometric phase is 



7r(uj + cjqA) 



7±,; 



: 7± - 7±,d 



IT 



<J r f A 



Based on Eqs. ([2]) and ([3]), we obtain a unitary gate with 
the dynamical and the Aharonov-Anandan phases with 
respect to the computational basis {|0), |1)}, 



V(r) 



cos 7 — 1 sin 7 cosx ~ ie ^ sin 7 sin x 
—ie 1 ^ sin 7 sinx cos 7 + i sin 7 cosx 
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FIG. 2: Schematic diagram of double- loop method for dynam- 
ical phase cancellation according to the proposal by Zhu and 
Wang 8] . Two time-dependent magnetic fields are applied se- 
quentially. The first magnetic field (loop 1) rotates counter- 
clockwise, while the second one (loop 2) rotates clockwise in 
order to eliminate a dynamical phase. 



where 7 = 
notation. 



2tt — 7+ = 7_ has been used to simplify the 



C. Cancellation of dynamical phases 

We closely follow Zhu and Wang's proposal [8[ in order 
to eliminate a dynamical phase. They proposed the use of 
two successive unitary operations, in which a dynamical 
phase cancels out while a geometric phase accumulates 
along these two operations. Each unitary operation asso- 
ciated with a loop is characterized by the time-dependent 
external field (|4]). The loop parameter corresponding to 
the ith loop is denoted by Qi(t) (i = 1, 2). Thus, in prin- 
ciple, we have four independent parameters in each loop, 
i.e., o^i, c^ r f, (pi, and uoi^. We note that they are not 
always tunable in a real experimental situation. 

We will search for the condition under which the dy- 
namical phases associated with the two loops cancel each 
other as shown in Fig.[2j We focus on the case in which 
(pi = and a;i 5r f = cc?2,rf(= o; r f) > for simplicity. The 
first loop (loop 1) is described by 

Fli(t) = (^1,1 coscjrft, — cc^i sinc<; r ft, u^o) , (7) 

while the second loop (loop 2) by 

^2(t) = — (^2,1 COSCJrft, — (J2,l SU1 LJ r f t, LV 2 ,o) Ry(@)- (8) 

Let R y (0) G SO (3) represent a rotation around the y- 
axis by an angle 0. The rotation angle 6 is chosen so 
that the corresponding cyclic vectors |^,±) for these two 
loops satisfies 




FIG. 3: (color online) Example of Zhu- Wang's double- loop 
method. The time-dependent external fields Qi(t) and Q2 (t) 
are shown in (a), while the closed trajectory on the Bloch 
sphere corresponding to the cyclic vectors 1^1,+) an d 1^2,+) 
in (b). We note that these are connected and thus form one 
closed trajectory. We set the loop parameters u;i,i = uj2,i = 
27T, o; r f = 0.7 x 27T, cji,o = 0.27 x 27T, and 002,0 = 1.5 x 2n in 
Eqs. (0 and J5J. We note that these parameters are calculated 
on the basis of a condition for nulling dynamical phases in 
Ref. [|[. In this example, T = \ in Eq. ([T3]) . 



where c G R. In other words, the Bloch vectors corre- 
sponding to \ipi,±) coincide with those to |^2,±)- Using 
the notation in Eq. ([5]), we find that in the loop i 

mi = kiOi, ki = (sinxi,0,cosxi), 

where 0\ = J3, O2 = Ry(0), and we write the 3x3 unit 
matrix as J3. The rotation angle \i is defined as 

tanxi = ^r^, Ai = cji,o - cj rf , (10) 
£±1 



tan% 2 



^2,1 
A 2 ' 



: ^2,0 + ^rf • 



(11) 



The angle is explicitly given as 
= X2~ Xi- 

Figure [3](a) shows an example of the time-dependent ex- 
ternal fields fti(t) and ^(t)- The corresponding closed 
trajectory on the Bloch sphere is drawn in Fig.GJb), in 
which the initial point corresponds to a cyclic vector 
l^i, +). Figure 02(b) shows that is n °t only cyclic 

for loop 1 but also for the total process (i.e., loop 1 and 
loop 2) due to the connection condition (|9j). It is neces- 
sary to search for oj^ a and uu r f (a = 0, 1) so that 



7i,d + 72,d = 0, 
7i,g + 72, g = Ttt. 

We write them more explicitly as follows: 



"2,1 



- CJ2,0 A2 



^2 



Ai A 2 

— - -\ = 2 • 



r, 



(12) 
(13) 



(14) 
(15) 



|^l,±>=^2,±>, 



(9) 



where Qi = {uxfi + A 2 ) 1 / 2 . There may be many com- 
binations of LOi^a and lo x { for a given T which satisfy the 
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Soft pulse 1 




Soft pulse 2 




Ry(9) 




Ry(-TT) 


Ry{-0) 



loop 1 



loop 2 



time 



FIG. 4: Schematic diagram of double-loop method for dynam- 
ical phase cancellation on the basis of a spin-echo approach. 
Two (four) soft (hard) square pulses are applied. We note 
that R y {6) = e - id(J y /2 , in which = X 2 - xi. 



conditions ([T4]) and ([T5j) [8]. We note that a set of the 
parameters employed in Fig. [3] is one example for the so- 
lution of Eqs. ([14]) and ([15]), in which r = ~. 

After the elimination of a dynamical phase, we have a 
one-qubit geometric quantum gate 



Vzw = 



AF-k 



iF-K 



l^i,-X^i 



(16) 



D. Spin-echo approach 



Zhu and Wang's proposal for eliminating a dynami- 
cal phase is not feasible for a conventional commercial 
NMR system where a field along the z-axis is strictly 
constant. In other words, it is difficult to realize H2 (t) 
in Eq. JSJ. In the present paper, we propose an exper- 
imentally feasible method, in which the loop 2 is di- 
vided into three successive steps while the loop 1 remains 
unchanged. The separation of the loop 2 is motivated 
by the spin-echo technique frequently employed in NMR 
experiments, in which the direction of time is reversed 
by an application of a pair of 7r-pulses. Three succes- 
sive operations are (a) a rotation around the y-axis by 
0(= X2 — Xi)> (b) an operation corresponding to preces- 
sion by a field — (^2,1 coscc; r ft, ^2,1 sinc^t, cj2,o) for a du- 



ration r = 27r/|cj r f |, and (c) a rotation around the ?/-axis 
by —0. The rotations R y (±0) correspond to the basis vec- 
tor change and back as shown in Fig.JU Rotation around 
the y-ax.is by is easy to realize by a radio- frequency (rf) 
pulse, which corresponds to the unitary operator 



Ry(6) 



9o y /2 



(17) 



We emphasize here that \ip2,±) = Ry(6)\ipi,±)- We 
assume that the pulse duration is infinitely short 
for simplicity. It should be noted that this opera- 
tion does not generate a dynamical phase since the 
?/-axis is perpendicular to both k\ and k 2 @, [Hj]. 
The operation corresponds to a precession by a field 
— (^2,1 coso; r ft, 002,1 sino; r ft, 002,0) for r = 27r/|o; r f| is 
given as 



u' 2 



in 2 T(-k 2 )-cr/2 



from Eq. J6j. The identity 



-k 2 )-(T/2 



implies that ~U' 2 can be realized by a precession under the 
field (c^2,i coso; r ft, — oj2,i sino; r ft, cc^.o) for r = 2it/\oj T f\ 
sandwiched by a pair of ±7r-pulses [25| . We again assume 
that R y (±.ir) is implemented for an infinitely short pulse 
for simplicity. 

Summarizing the above arguments, the total process is 
described by R y {-0)R y {-TT)U 2 {T)Ry{TT)R y {0)U 1 {r) 0, 



which is equivalent to Eq. ([16]) if the dynamical phase is 
zero. The unitary operator Ui(r) is defined as Ui(r) = 
e iu T iT(jy/2 e in i ki-cT/2 _ rp^ e g eome ^ c g a ^ e w hich we are 

going to demonstrate takes the form 



echo 



Ryi-Tr^^Ryi^RyiOW^T). (18) 



Although the complete realization of Zhu and Wang's 
original proposal ([T6]) requires R y (—0) at the end of the 



process ([T8]h we can omit it for constructing a geometric 
quantum gate since R y (—0) does not generate any dy- 
namical phase here. We note that Vzw = Ry(— 6)U ec ho- 
Under the conditions ([T2]) and ([T3]h the matrix represen- 
tation of Eq. (jTHJ) in the computational basis {|0), |1)} is 
given by 



cos(Ftt) cos((9/2) 
cos(rTr) sin(<9/2) 



i sin(r7r) cos 6 — cos(r7r) sin(0/2) — i sin(r7r) sin 6 
i sin(r7r) sin O cos(r7r) cos(#/2) + i sin(r7r) cos O 



(19) 



r 



where = (x2 + Xi)/2- We note that r-zr is the total 
geometric phase. The pulse sequence JTHJ leads to intu- 
itive understanding of the cancellation mechanism of the 
dynamical phase. Let us consider the case of xi — X2> 



i.e., the two loop are completely identical, for simplicity. 
If no 7r-pulse is applied, the dynamical property of the 
loop 1 is the same as that of the loop 2 and the total dy- 
namical phase is the addition between the contributions 
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from the loops 1 and 2. It should be noted here that the 
±7r-pulses induce the time-reversal dynamics in the form 
of U' 2 in the loop 2. Under the time-reversal transfor- 
mation, the sign of the dynamical phase associated with 
the loop 2 is inverted and hence the dynamical phase is 
completely eliminated. It is necessary to employ differ- 
ent processes between the loops 1 and 2 to prevent the 
cancellation of the geometric phase associated with the 
two loops. The matrix representation (fT9|) implies that 
C4cho contains three parameters T, and O. Due to the 
limitation in the control parameters, it may be difficult 
to choose them independently in a standard liquid-state 
NMR. We will show that 9 and should be regarded 
as functions of T and uoi in Sec. Ill E[ in order to satisfy 
Eqs. ([12]) and (13]) as shown in Table H On the other 
hand, we are still able to use the rf phase <j). Thus, we 
have the necessary number of free parameters to express 
arbitrary elements of SU(2). 

E. Implementation in liquid-state NMR 

We implement the double-loop scheme in liquid-state 
NMR. We take different loop parameterization from that 
of Zhu and Wang for ease of implementation. We con- 
sider the system in a rotating frame defined by o; re f. 
Hereafter, we will denote an angular frequency x mea- 
sured with respect to this rotating frame as x = x — o; re f . 
Thus, one can explicitly understand which quantities are 
controllable by choosing a proper rotating frame. 

We take a common value 

U = (JO - ^ref < 

to coi^o and 002,0- The value of ujq in the experiment 
will be explained in Sec. lIIIBl Similarly, we assume that 
^1,1 = wi and a;2,i = Instead of these simplification, 
we allow different values with respect to uo T { between the 
two loops, i.e. cj r f = a;i >r f in Eq. and u r f = cj2, r f in 
Eq. J8|). These changes do not alter the basic strategy 
for searching parameters that satisfy 71 ^ + 72, d = and 
7i,g + 72,g = T7T. We consider the two loops in the rotat- 
ing frame in which the frequency is a)^ r f, the amplitude 
(Ji, and the duration 7^ = 27r/|o>^ jr f | , i.e. 

£ti(t) = (tjJl COSCc^rf t, — (J 1 Sinc^rf t, UJq) (0 < t < Ti) . 

The solutions ^, r f/|^o| are numerically obtained for 
given e(= ooi/\ujo\) and T. 

We summarize our parameter choice. First of all, we 
adopt a common value to uo\^ and 6^2,0, i-e., Cjo in the 
rotating frame defined as the angular frequency (j re f . The 
value of uj\{= = (^1,2) is given by uo\ = e|tDo|? in 
which e is a positive number. For a given e and an aimed 
geometric phase T, we can numerically find proper a)^ r f so 
that 7i 5 d + 72, d = and 7i >g + 72, g = Ttt. The results for 
r = I and e = 0.5, 0.3, and 0.1, for example, are shown 
in Table [H From the observation of Eqs. ([TO]) and (fTT]h 
the sign of u)i jr f should be opposite to the one of 0)2, rf- 



It should be noted that the parameters given in Table 
|T] are compatible with this requirement. The resultant 
geometric quantum gate is Eq. ([T9]h The values of 6{= 
X2 — Xi) an d ©(= (X2 + Xi)/2) are given in Table H 
When r = |. Eq. (fT9|) takes the form 

^echo(O) = e- i7r ^ /2 e^ 2r2fc2,<T/2 

Xe -i(e+7r)a y /2 e in 1 T 1 k 1 -cr/2 

= e _ W2 /cos9 sinO\ 

y sin 6 — cos By' v 7 

which we experimentally demonstrate in the next section. 

III. EXPERIMENTS 
A. Sample and spectrometer 

We implement a one-qubit gate described by Eq. ([20]) 
with a conventional commercial NMR system. We em- 
ployed a JEOL ECA-500 NMR spectrometer whose 
hydrogen Larmor frequency is approximately 500 MHz. 
13 C nucleus in a 0.6 ml, 0.2 M sample of 13 C-labeled chlo- 
roform (Cambridge Isotope) in d-6 acetone is employed 
as a qubit, while protons are decoupled by a standard de- 
coupling technique, called WALTZ [25]. We have chosen 
13 C-labeled chloroform for future experiments involving 
two-qubit gates. The transverse and the longitudinal re- 
laxation times are T2 ~ 0.3 s and T\ ~ 5 s, respectively. 
The longitudinal relaxation time is shorten by adding a 
small amount of Iron(HI)-acetylacetonate so that a repe- 
tition rate can be increased. T2 and T\ without Iron(IH)- 
acetylacetonate are ~ 0.3 s and ~ 20 s, respectively. 

B. Pulse sequence 

As we discussed in the previous section, the gate ([20]) 
can be realized with two rotating magnetic fields and two 
hard (short) pulses. The rotating fields are effectively ob- 
tained by two soft (long) pulses which are rotating with 
different frequencies o)^ r f = cc^ r f — co re f (i = 1,2) in the 
rotating frame with frequency cj re f. The first soft pulse 
(loop 1) is a usual square pulse, while the second soft 
pulse (loop 2) is a (frequency) shifted laminar square 
pulse (SLP) [31] . This SLP is employed in order to ob- 
tain the same phase <p in Eq. ^ for loop 2 as that for 
loop 1, i.e., 0i = 2 . 



e 


£l,rf/&0 


^2,rf/^0 


[rad] 


6 [rad] 


0.5 


-0.6815 


0.7803 


-0.7298 


-0.639 


0.3 


-0.8221 


1.105 


-0.9571 


-0.589 


0.1 


-0.9422 


1.609 


-1.008 


-0.542 



TABLE I: Solutions of Eqs. (T2J and (H3l) for T = \ and e = 
^i/po| = 0.5, 0.3, 0.1 in the rotating frame so that (Do < 0. 
We note that Ti\Cji fT f | = 27T, 6 = X2 — Xi> an d B = (xi+X2)/2. 
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e 


F e (Io,£oW-£ ) 


Tr[S(I )] 


F e (I ,£ 2 ) 


Tr[S^Io)] 


0.5 


0.75 


1.00 


0.74 


1.02 


0.3 


0.88 


1.08 


0.83 


1.07 


0.1 


0.84 


1.07 


0.85 


1.06 



TABLE II: The entanglement fidelities for single and double 
operations with e = 0.5, 0.3, and 0.1. 

We take \uj \ = 2tt x lOOOrad/s and fa = through- 
out the experiments. The condition fa = is taken for 
simplicity as mentioned in the beginning of Sec. Ill CI We 
independently calibrate the strengths of the soft and hard 
pulses in order to minimize a non-linearity error in set- 
ting the rf pulse amplitude. The duration of a hard 
7r-pulse is set to 21.6 /is throughout the experiments. We 
ignore t^p in setting the phase of the second soft pulse, 
which is justified by the fact that thppol <C 2tt. The pre- 
cision of pulse duration control is 100 ns. The durations 
U ySp of two soft pulses are set to 

U,sp\^i,ri\ = 27T. 

We demonstrate three different gates with e = 0.5,0.3 
and 0.1. We note that the phase of the second hard pulse 
corresponding to R y (—7r) must be adjusted, presumably 
because the oscillator in the NMR spectrometer is dis- 
turbed in generating a SLP. It should be recalled that a 
SLP employs intensive phase modulation. 




C. Results 

Implemented gates with T = | are evaluated by per- 
forming quantum process tomography [32] . The practical 
details are explained in Ref. [331 ]. A quantum process £, 
such as a gate operation or relaxation process, is 



Efc E l E k 



in the operator sum (or Kraus) representation [34], |35[ . 
When all E^s are determined, £ is considered to be iden- 
tified. This identification is called quantum process to- 
mography. 

The Bloch sphere in Fig.[5](a) is mapped under the gate 
operations to the surfaces in Figs.EJb), (c) and (d), which 
correspond to e = 0.5,0.3 and 0.1, respectively. If the 
gate operations are perfect, the surfaces are the spheres 
of unit radius (i.e., the Bloch sphere). The left panel of 
each row shows the theoretical final state. In the middle 
panels, the results for the single gate operation are shown. 
Finally, the right panels are for the two-successive (dou- 
ble) gate operation. The Hadamard gate obtained when 
O = — 7r/4 is, for comparison, shown in the right pannel 
of Fig.[5ja). From these figures, we find that /7 e cho(@) in 
Eq. ([20]) is implemented although it is not perfect. 

We numerically evaluated the fidelity of the imple- 
mented gate using the entanglement fidelity [H, [35[ given 




FIG. 5: Gate operations visualized. The Bloch sphere in 
(a) is mapped to the surfaces in (b), (c), and (d) under the 
gates with e = 0.5, 0.3 and 0.1, respectively. The right surface 
in (a) is an expected Bloch sphere when G = — 7r/4, which 
corresponds to the Hadamard gate. Each left panel in (b), 
(c), and (d) corresponds to the theoretical final state. The 
middle panels are the results for the single gate operation 
£4cho(®)- The right panels are for the two- successive (double) 
gate operations. 



by 

F(I Sou-' ) = E fc |Tr[E fc £/ echo (-9)Jo]| 2 

e( °' echJ Tr[E k E k u ec a-e)ioUl ho (-e)Et] 

= E fc |Tr[£ fc £/ echo (-9)I ]| 2 
Tr(J2 k E k I El) ' 

where U ec ho is a super operator corresponding to 
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the unitary operator U echo (Q) (i.e., U echo (p) = 
^echo(@)p^ e t cho(®))' ^0 = 1/2, and 1 is the identity ma- 
trix of dimension 2. One can find that F e (/o,£o^ ) = 
1 when the gate operation S is perfect. In the case of 
two successive gate operation, F e (io, £ 2 ) gives a measure 
of the fidelity since [/7 ec ho(@)] 2 = — 1. The entanglement 
fidelities corresponding to the gate operations are sum- 
marized in Table HT1 

The fidelities of the demonstrated gates are not high. 
This may be attributed to the inhomogeneous rf field. 
The free induction decay signal of the thermal state af- 
ter a 57r/2-pulse, which corresponds to the operation 
e -57va x /4 £ Qr exam pi e5 reduces to about 85% of that after 
a 7r/2-pulse, which corresponds to the operation e -7 ^/ 4 
for example. This fact indicates that there is some rf 
field inhomogeneity which may account for most of the 
reduction in the fidelities in Table HI1 Pulse sequences in 
usual NMR operations are designed so that the rf field 
inhomogeneity does not affect measurements, for exam- 
ple, by employing composite pulses. Such techniques are 
not available in our experiments. 

IV. SUMMARY 

We demonstrated the elimination of the dynamical 
phase and the implementation of the quantum gates 



with pure nonadiabatic geometric phases in a liquid- 
state NMR quantum computer, based on the double-loop 
method. By means of a spin echo technique, we modified 
the original proposal so that quantum gates are imple- 
mented in a standard high precession NMR system for 
chemical analysis. We have proposed and experimentally 
verified an alternative method to eliminate dynamical 
phase. The extension of the present method to two-qubit 
operations is an important future work [36[ . We believe 
that our work is the first step toward physical realization 
of working geometric quantum gates and further efforts 
should be made for improvement of the gates. 
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